Scaling of magnetic fluctuations near a quantum phase transition 
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We use inelastic neutron scattering to measure the magnetic fluctuations in a single crystal of 
the heavy fermion alloy CeCus.gAuo.i close to the antiferromagnetic quantum critical point. The 
energy(E)-, wavevector(q)- and temperature(T)-dependent spectra obey E/T scaling at Q near 
(1,0,0). The neutron data and earlier bulk susceptibility are consistent with the form ~ f(Q) + 
(— iE + aT) a , with an anomalous exponent a ~ 0.8 7^ 1. We confirm the earlier observation of 
quasi-low dimensionality and show how both the magnetic fluctuations and the thermodynamics 
can be understood in terms of a quantum Lifshitz point. 

PACS numbers: 71.27.+a, 75.20.Hr, 75.40.Gb 



Quantum phase transitions are zero temperature tran- 
sitions driven by a parameter such as pressure, mag- 
netic field, or composition which regulates the ampli- 
tude of quantum fluctuations. An important realization 
of the last years has been that many of the most in- 
teresting phenomena in condensed matter physics seem 
to occur near such transitions. Notable n examples in- 
clude high-temperature superconductivitycra and heavy 
fermion behavior.EHa For high-temperature (cuprate) su- 
perconductivity, the quantum critical hypothesis leads 
very naturally to one of the key attributes of the normal 
state, namely that the energy scale governing spin and 
charge fluctuations isJjhe temperature T itself, a property 
labelled E/T scaling^ 

Heavy fermion behavior is generally understood as a 
consequence of a competition between intersite spin cou- 
plings and the single-site Kondo effect Antiferromag- 
netic order develops in a heavy fermion metal at a quan- 
tum critical point where the intersite spin couplings over- 
come the disordering Kondo effect of local spin-exchange 
processes. While the notion of such competing inter- 
actions is old, the nature of the quantum phase tran- 
sition is a subject of great current interest. CeCug, a 
heavy fermion_metal with one of the largest known linear 
specific heatsJ3 provides an ideal opportunity to explore 
this physics in quantitative detail. Here, the quantum 
critical point (QCP) occurs as Au is substituted for the 
Cu atoms; when more than 0.1 Cu sites per Ce are re- 
placed, the heavy fermion paramagnet gives way to an or- 
dered antiferromagnet.il Fig. [I] illustrates the associated 
phase diagram and ordering vectors. In spite of the lagse 
and rapidly growing literature on this particular QCPtffl, 
the E/T scaling which has been so thoroughly tested 
for the single crystal cupratestl and for polycrystalline 
UCus-^Pd^EI has not yet been observed for single crystal 
heavy fermion systems in general and CeCug-zAua; in 
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FIG. 1. Location of the magnetic Bragg peaks 

or spin fluctuation maxima r^Q^i CeCue-zAu^ at different 
Au-concentration x (from refsLTtJ). Dotted lines are the tra- 
jectories in the (hkO)-plane of the constant-E scans presented 
in Fig. EL Inset shows magnetic ordering temperature Tjv 
(from rem) vs. x. The arrow marks the thermal trajectory 
investigated for x=0.1 in this paper. 

particular. Our purpose here is to report the discovery, 
obtained from magnetic neutron scattering, of E/T scal- 
ing near the QCP in CeCus.gAuo.i. Beyond establishing 
the universality of E/T scaling among QCPs hosted by 
vastly different alloys, our data yield a non-mean-field 
value for a key critical exponent characterizing the quan- 
tum phase transition in CeCus.gAuo.i. 

We used the Czochralski technique to grow a single 
crystal of CeCus.gAuo.i with volume ~4 cm 3 and a mo- 
saic of 2.5°. To establish the homogeneity and com- 
position of the sample, we used high-resolution neutron 
diffraction in the orthorhombic [001] zone employed for 
our inelastic neutron scattering measurements. The re- 
sult was that throughout the whole sample, the room 
temperature (Pnma) lattice constants were a = 8.118(2) 
A and b = 5.100(1) A corresponding to the Au con- 
centration of x= 0.10±0.03.B Furthermore, the transition 
into the low-temperature monoclinic phase occurred at a 
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FIG. 2. Constant-E scans through the magnetic scattering 
from CeCus.gAuo.i at Eo=0.167 meV near Q=(0. 8,0,0) (a,b) 
and Q=(0 1±0.15 0) (c,d). Solid lines in (a,b) correspond 
to Eq.(3), corrected for experimental resolution and with the 
same parameters as shown in Fig. | and f(Q) expanded in 
even powers in Q. The lines in (c,d) are simply guides to the 
eye. The T-independent peak at (010) in (c) is of nuclear 
origin. 



location- independent T5 = 70±2 K, again confirming the 
homogeneity of the sample. Between T5 and 5K, the an- 
gle 7 grows from 90° to 90.7°, while the a-axis direction 
remains unchanged. Thus, the monoclinic distortion is 
sufficiently small that we follow past custom and use or- 
thorhombic notation to label points in reciprocal space. 

After installing the sample in the Bell Laborato- 
ries/Ris0 dilution refrigerator, we collected inelastic neu- 
tron scattering data using TAS7 at the Ris0 DR3 reactor 
equipped with a PG(002) monochromator and analyzer 
and a BeO filter. The final neutron energy was fixed 
at E/= 3.7 meV. The energy resolution, defined as the 
full-width-at-half-maximum (FWHM) of the elastic inco- 
herent scattering from the sample, was 0.134 meV. 

When we began our experiments, we knew that the 
most intense magnetic fluctuations in the paramag- 
netic parent CeCug were located close to the inoquiv- 
alent reciprocal lattice points (010) and (100). □ At 
the same time, the incommensurate vectors (l±5 0) 
fully describe the magnetic order in CeGue-xAurr with 
x>0.5BB and do so partially for x= 0.2.113 Thus, to dis- 
cover the dominant fluctuations near the QCP at x= 0.1, 
we investigated fluctuations with wave-vectors of both 
types. Fig. shows the corresponding constant-energy 
scans at Eo= 0.167meV, collected along the trajectories 
indicated in Fig. |l|. Frames (c) and (d) demonstrate that 
incommensurate peaks survive at (0 1±0.15 0), and that 
their intensity grows upon cooling. The widths of the (0 
1±0.15 0) peaks in the (1,0,0) and (0,1,0) directions are 
similar and not resolution-limited, corresponding to mag- 
netic correlations with a range of order 8A(at this Eo). 
Even though they are still present and growing with de- 
creasing T, the (0 1±0. 15 0) peaks are five times weaker 
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FIG. 3. Constant-Q scans at Q=(0.8 0) and Q=(1.8 0) 
and (0.6 0.3 0) at different temperatures T vs energy tranfer 
E (a) and at different E vs T (b). Lines correspond to Eq(3) 
with c=83cts/8min-meV ' 74 , where a and a are given by the 
analysis of the scaled data (see Fig. ^ and text) and include 
resolution corrections. The optimal (T- and E-independent) 
values for f(Q) are 0.036(4)meV ' 74 and 0.38(2)meV 74 for 
Q=(0.8 0) and the other Q-values respectively. Note that 
f(0.8 0) is finite, indicating that our sample is slightly sub- 
critical. But it is sufficiently small that when the data are 
examined outside the regime where the resolution matters, 
the deviations from criticality are negligible. 



than the scattering, shown in frames (a) and (b), near 
(100). Here, there is a weak modulation along (1,0,0) 
with hints of maxima at (1±S 0) for S » 0.2, and a 
much sharper modulation along (0,1,0). Both the peak 
intensity and width along (0,1,0) are clearly T and E- 
dependent, corresponding to growing antiferromagnetic 
correlations with decreasing T and E. Furthermore, T 
and E appear interchangeable in the sense that the pro- 
files for E=0.3meV > k B T=0.006meV and E=0.167meV 
< ksT=0.345meV are indistinguishable. That the larger 
of E and k^T controls the width k(E,T) is expected for 
quantum critical points and has been previously observed 
in the context of high-temperature superconductors .□ 

We next characterize the spectrum of the most intense 
fluctuations, namely those at Q=(l±0.2 0). Fig. || 
shows such constant-Q spectra for several T. For com- 
parison, it also shows data at Q=(1.8 0) and (0.6 0.3 
0), for which no special Q- or T-dependent enhancement 
is found. The Q- and T-independent elastic scattering 
background has been subtracted from all data shown. At 
Q=(0.8 0), the spectra narrow dramatically with reduc- 
tion in T, to the point where they have a rising edge at 
low E which is comparable to that of an instrumentally 
broadened step function. Thus, as happens near any sec- 
ond order phase transition, cooling results in a reduced 
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magnetic relaxation rate. Even so, a long high-E tail of 
the spectra persists at all T that the FWHM never ap- 
proaches that of the instrumental resolution. At Q=(1.8 
0), the spectra change much less dramatically. 

In addition to performing measurements as a function 
of E at fixed T, we have also collected data as a function 
of T for fixed E (Fig. |J). The similarity between the 
T-scans (b) and the E-scans (a) suggests E/T scaling, 
which we have tested by fitting the data to the form 

X "(E,T) = T- a g(E/k B T). (1) 

No assumptions were made about g(x) beyond the ability 
to approximate it by a histogram with stepsize 0.1 on a 
logio E/T scale between -0.5 and 0.7. The best collapse 
of the data onto a single curve, as checked by the small- 
est (log) deviation (<Jiog(g) see inset) from the mean step 
values, is obtained for a = 0.75 ± 0.1 and shown in Fig. 
|], which displays all the data at Q=(0.8 0) unaffected 
by the finite experimental resolution. 

Having confirmed E/T scaling, we ask whether the 
associated exponent a is consistent with other data and 
theory. The simplest mean field approachli3 dictates a 
susceptibility \ 

X {Q,E,T)- l = c- 1 U{Q)-iE + aT), (2) 

where f(Q) is a smooth function of Q with minima at the 
wave- vectors characterizing the eventual magnetic order. 
If we position ourselves at the critical Q and composition, 
f(Q)=0 and x" — c E/[(aT) 2 +E 2 ], which we can compare 
to Eq.(l) to read off the scaling function g (y) = y/(l + 
y 2 ) with y = E/aT and the exponent a=l. The dotted 
line in Fig. ||(a) corresponds to the (log) optimized g (y) 
with a/ke=l and deviates clearly from the data. 

Neither g (y) nor a= 1 account for the measurements; 
thus the simple mean-field theory fails for CeCus.gAuo.i. 
To make progress, we modify Eq(2) to account for arbi- 
trary values of a via the simple expedient of replacing 
(-iE + aT) in Eq.(2) by (-iE + aT) a , 

X- 1 =c- 1 {f{Q) + (-iE + aT) a ) (3) 

Beyond accounting for the anomalous value of the expo- 
nent of T in the scaling relation, this generalization has 
two other important implications. The first is that the 
scaling function (for Q at which f(Q)= 0) becomes 

g(y) = csinlatan- 1 ^)]/^ 2 + l) a ' 2 (4) 

The solid line in Fig. ^ corresponds to the best fit of 
Eq.(4) to the data (with c=88cts/8min-meV - 74 ); note 
that the value obtained from floating a in Eq.(4) in the 
fitting process is at 0.74±0.1 indistinguishable from that 
obtained as the exponent for the thermal prefactor a in 
Eq.(l). At the same time, a/fcs=0.82 is a number close to 
unity, indicating that the magnetic fluctuations have an 
underlying energy scale very close to the thermal energy. 

The second consequence of the generalized form (3) is 
that for E=0, x'iQ^)' 1 = ^(/(Q) + (aT) a ), which 
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FIG. 4. (a) Scaling plot of the data from Fig. | for Q=(0.8 
0). The inset shows how the quality of the scaling collapse 
varies with a in Eq.(3), using a procedure unbiased towards 
any particular form of the scaling function g(y) (see text). 
Solid and dashed lines correspond to Eq.(4) with a=0.74 and 
1, respectively, (b) Double logarithmic plot of susceptibility 
Xbuik (from M/B in B c =0.1T froml) vs. T after subtract- 
ing a zero T value xo corresponding to a 'parallel' shunt, as 
suggested by Eq.(3). 

implies that {x'iQXt 1 - x'{Q, 0)" 1 ) -1 = c(aT)-«.We 
test this relation by using the very accurate and dense 
data provided by bulk magnetometryu for a special value 
of Q, namely Q=0. The fact that on the double loga- 
rithmic scale of Fig. |J(b) the data lie on a straight line 
shows that the relation is indeed satisfied. The line is 
not parallel to the diagonal for which a=l, but indeed 
has slope a=0.8±0.1, indistinguishable from a obtained 
from the neutron scattering data by the other two meth- 
ods described above. 

The neutron results are consistent with the bulk sus- 
ceptibility: do they also account for the thermodynam- 
ics? One of the key features of this system, is that, at low 
T, C/T = 7(T) diverges logarithmically.B To see whether 
j(T) ~ ln(T) is consistent with the neutron measure- 
ments, we need to consider the detailed dependence of 
f(Q). In keeping with earlier investigationsQ, the criti- 
cal fluctuations arc almost two-dimensional. Specifically, 
the constant E scans along (1,0,0) (Fig. |^(b)) look flat- 
topped, indicating that the quadratic term in the power 
series expansion of f(Q) around (100) almost vanishes in 
this direction. A good fit to the data in the vicinity of 
Q = (0.800) is given by f((100)+SQ) = (0.036+12(5/c) 2 - 
0.08((5/i) 2 +2((5/i) 4 )meV - 74 (see dashed/solid lines in Fig. 
||). That the quadratic term Bh(Sh) 2 is negative (with 
a positive quartic term Ch(5h) 4 ) also follows from the 
observation that CeCue-xAua, orders not at (100) but 
at pairs of incommensurate wavevectors nearby (see Fig. 
[l]). As the sign of B/j changes, the minimum at 8Q = 
of f(SQ) splits into the pair SQ = ± v / -B h /2C h . In 
the classical theory of phase transitions, vanishing of the 
quadratic stiffness occurs at a Lifshitz pointta. Here, we 
are dealing with a quantum Lifshitz point. 

To understand the thermodynamics of a quantum Lif- 
shitz point, we assume that the quadratic stiffness van- 
ishes along only one direction (e.g. (1,0,0)) in reciprocal 
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space or more generally, that it vanishes for a finite num- 
ber qf-|lines centered at certain points Q in reciprocal 
spaceu-3. The most important contributions to \ will then 
be proportional to 1/[B ± (SQ ± ) 2 + C||(<5Q||) 4 + {-iE + 

aT)~] in the vicinity of Q, where 8Qu and <5Qj_ measure 
the departure from Q parallel and perpendicular to the 
line direction, and we have employed the notation of dy- 
namical critical phenomena, writing a — 2/z. If k± and 
Ky are the inverse correlation lengths in these directions 

so the dynamical critical exponent is 2z along the (1, 0, 0) 
direction, but z in the perpendicular directions. This 
has the effect of reducing the effective dimensionality 
D e //. For an isotropic QCP, the free energy scales as 
F ~ Tk d ~ yi+^/z p or a q uan tum Lifshitz point, we 
must now write 

F ~ T/t^ _1 K|| ~ 7 ll +(- D_ 5)/ z . 

So D e f f is lower by one half, and 

7 (T) = -d%F(T) ~ T< D -iW2-i. 

We see that for D = 3, a = 0.8, the exponent for j(T) 
vanishes, in accord with the specific heat data. If we 
ignore the anisotropic form of f(Q) indicated by the neu- 
tron scans through reciprocal space and assume a con- 
ventional (non-Lifshitz) QCP for CeC^-zAu^, the a ob- 
tained from the energy- and temperature-dependence of 
the neutron data as well as the bulk \ would lead to 
j ~ T°- 2 and T~ 0/2 for D=3 and 2 respectively, in dis- 
agreement with the C(T) results. 

Our ability to use the phenomenological scaling form 
for x _1 for both the thermodynamics and the uniform 
susceptibility indicates that the anomalous energy depen- 
dence of the susceptibility extends over a wide region of 
the Brillouin zone. We can write the T=0 susceptibility 
in (3) as x~HQ,E) = X^{E) - J(Q), where X {E) de- 
scribes the response of individual local moments and sur- 
rounding conduction electrons and J{Q) defines the inter- 
actions between different screened local moments. Eq.(3) 
actually determines \^ 1 {E) and J(Q) up to the same ad- 
ditive constant O which sets the scale for the screening: 
x -!(£) = c - x {-iE) a + 9 and J(Q) = - C -\f(Q) + 9. 
Our result that a = 0.8 ± 0.1 implies that x(-E') is non- 
analytic near E=0. In the "Millis-Hcrtz" theory of a 
quantum critical point, the soft magnetic fluctuations 
couple to a Fermi liquid and the corresponding local re- 
sponse is analytic, with a = 1. The anomalous exponent 
a / 1 for CeCus.gAuo.i thus represents a fundamental 
iocai deviation from Fermi liquid behavior. One possi- 
ble interpretation is that the magnetic QCP has quali- 
tatively modified the underlying moment compensation 
mechanism,t3 causing the local electron—propagator to 
develop an anomalous scaling dimension £3 

We have performed the first detailed single crystal 
measurements of the magnetic fluctuation spectrum in a 



heavy fermion system close to an antiferromagnetic quan- 
tum critical point. We find that the data obey E/T scal- 
ing with an anomalous scaling exponent a « 0.8. The 
same exponent appears in three other independent quan- 
tities: the scaling Junction itself, previously published 
bulk susceptibility|j and the specific heat. To obtain 
the scaling relation connecting C(T) with the magnetic 
measurement, we take into account the (again indepen- 
dent) indications of our neutron data that the QCP in 
CeCus.gAuo.i is actually near a quantum Lifshitz point, 
whose hallmark is the near degeneracy of ordered states 
with characteristic wavevectors along some line(s) in re- 
ciprocal space. The general idea that there are many 
potential ground states for CeCus.gAuo.i is further rein- 
forced by our observation of an enhanced susceptibility 
at the seemingly inequivalent locations (0 1±0.15 0) and 
(1±0.2 0) in reciprocal space. 
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